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Abstract 

A reaction-diffusion problem with an obstacle potential is considered in a bounded 
domain of R^. Under the assumption that the obstacle ^ is a closed convex and bounded 
subset of R" with smooth boundary or it is a closed n-dimensional simplex, we prove 
that the long-time behavior of the solution semigroup associated with this problem can be 
described in terms of an exponential attractor. In particular, the latter means that the 
fractal dimension of the associated global attractor is also finite. 



1 Introduction 

This paper is devoted to the long-time behavior of solutions of the foUowing reaction-diffusion 
system with an obstacle potential in a bounded and regular domain fl C 

(dtu- A^u + dljt{u) - Xu3 0, ^^^^ 

Here u = x), ■ ■ ■ ,Un{t,x)) is an unknown vector- valued function, A^. is a Laplacian with 

respect to the variable a;, A > is a given constant, .J(^ is a given bounded closed convex set in 
R" containing zero and dlje stands for the subdifferential of its indicator function Ij^: 

T / N 1 0, u e ,Jf, , 

^^(") := , ^ (1-2) 

I oo, u f . 

We recall that the subdifferential dlj^ consists in the set of vectors in K" such that y G dlj(^{x) 
if and only if (y, x — w)^^ > for any w G J^, where (•, •)rii is the scalar product in R". It is 
also well known that dlje turns out to be a multivalued maximal monotone operator (see [6] , 
pg.25; see also Section [5] below for the rigorous definitions). 

Equations and systems of the type (jl.ip appear quite often in the mathematical analysis 
of phase transitions models or in reaction diffusion processes with constraints. In the first 
physical situation, (jl.ip rules the evolution of the so called order parameter, which is an n- 
dimensional vector in the case of multicomponent systems (see, e.g., [7]). Moreover, since the 
order parameter u is usually related to the pointwise proportions of the n independent phases 
shown by system under study, it is physically reasonable the fact that it attains values only in 



a bounded (convex) subset of M", usually an n-dimcnsional simplex 

JT := |p = (pi, ... ,p„) e R" such that < 1, > 0, i = 1, • • • , n| . (1.3) 

This is motivated by the requirement that no void nor overlapping should appear between the 
phases. In particular, equation (jl.ip with n = 2 and A = appears in the Fremond models of 
shape memory alloys (see [El also for other models of phase change showing the ubiquity of 
subdifferential operators in this framework) where the two components ui and U2 denote the 
pointwise proportions of the two martensitic variants. Finally, in the simpler scalar case, i.e. 
n = 1 and J(f ~ [0, 1] (for instance), equation (jl.ip is usually referred as Allen Cahn equation 
with double obstacle. 

The mathematical analysis of equations of the type (jl.ip (and more general equations 
associated with maximal monotone operators in Hilbcrt spaces) has attracted the attention of 
researchers for many years. In the particular case of equation p.ip . results concerning existence, 
approximation and long time behavior of solutions (e.g., in terms of global attractors) are known 
and by now classic (without any sake of completeness and referring only to results in the Hilbert 
space framework, we quote [S], [5], [3], [IE], [^). 

However, to the best of our knowledge, the finite/infinite dimensionality of the global 
attractor associated with the obstacle problems has not been yet understood (even in the 
simplest case of AUen-Cahn equation with double obstacle). Indeed, the classical machinery for 
proving the finite-dimensionality (in terms of fractal or/and Hausdorff dimension) of the global 
attractor (which perfectly works in many cases of dissipative systems generated by non-linear 
PDEs with regular non-linearities, see [21 [28] and references therein) is based on the so-called 
volume contraction arguments and requires the associated solution semigroup to be (uniformly 
quasi-) differentiable with respect to the initial data at least on the attractor. 

Unfortunately, this differentiability condition is usually violated if the underlying PDE 
has singularities or/and degenerations and, in particular, it is clearly violated for the obstacle 
problems like ()l.ip . Thus, the classical scheme is not applicable here and this makes the problem 
much more difficult and interesting. In fact, it has been recently shown that, contrary to the 
usual regular case, the singular/degenerate dissipative systems can easily generate infinite- 
dimensional attractors even in bounded domains. For instance, the global attractor of the 
degenerate analogue of the real Ginzburg-Landau equation 

dtu = A,{u^) + u - u\ u\^^=0 (1.4) 

is infinite-dimensional for any bounded domain of (thanks to the degeneration at u = 0), 
see [TT] . On the other hand, in recent years the finite-dimensionality of the global attractor result 
has been established for many important classes of degenerate/singular dissipative systems 
including Cahn-Hilliard equations with logarithmic potentials (see [22)). porous media equations 
(under some natural restrictions which exclude the example of (|1.4p . see |llj). doubly non- 
linear parabolic equations of different types (see [SH [25l [12] ) , etc. In these papers, the finite 
dimensionality of the global attractor is typically a consequence of the existence of a more 
refined object called exponential attractor, whose existence proof is often based on proper forms 
of the so called squeezing/smoothing property for the differences of solutions. 

We remind that the concept of exponential attractor has been introduced in [S] in order 
to overcome two major drawbacks of the global attractors: the slow (uncontrollable) rate of 
attraction and the sensitivity to perturbations. Roughly speaking, an exponential attractor 
(which always contains the global one) is a compact finite- dimensional set in the phase space 
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which attracts exponentially fast the images of all bounded sets as time tends to infinity (see 
Section [3] for the rigorous definition). Thus it turns out that, in contrast to global attractors, the 
exponential attractors are much more robust to perturbations (usually, Holder continuous with 
respect to the perturbation parameter). Moreover, the rate of convergence to the exponential 
attractor can be controlled in term of physical parameters of the system considered, see |H] and 
the more recent survey |23| (and also references therein) for more details. Finally, the finite 
dimensionality of the global attractor immediately follows from the finite dimensionality of the 
exponential attractor. 

The main aim of the present paper is to extend the exponential attractors theory to some 
classes of reaction-diffusion problems with obstacle potentials. Although the methods based on 
the proper squeezing/smoothing property for the differences of solutions do not require the 
differentiability with respect to the initial data and, in principle, can be applied also to the 
obstacle problem (jl.ip . this application is far from being straightforward in our situation since, 
to this end, one needs to produce estimates for the difference between the Lagrange multipliers 
(namely the selection of the subdiffcrcntial dlj^ (u) which turns the differential inclusion p.ip 
into an equation) associated with two solutions ui{t) and U2{t). This kind of estimates, which 
form the Assumption ^ (sec p.20p and p.lOp for the rigorous formulation), roughly speaking 
look as follows 

/ \\dly(r{ui{t))~dly^{u2m\LHn)dt<C\\ui{0)-U2m\L^ (1.5) 
Jo 

where ui and U2 denote two solutions starting from the proper absorbing set and, with a little 
abuse of notation, we refer to dijfr as it were single valued. Such kind of estimates, to the best 
of our knowledge, do not seem to be already known. 

This paper is organized as follows In Section [2 we recall the basic results related with the 
wellposedness and regularity of solutions of the obstacle problem As usual, the results 

on the singular problem (|l.ip are obtained by approximating the singular obstacle potential 
by more regular ones. The usual regularization used for this kind of problems is the Moreau- 
Yosida approximation (see, e.g., [B]). However, it turns out that, in order to prove estimates 
of the type (jl.Sp . it is more convenient to implement different kind of approximation schemes. 
Thus, in Section [2] we also give a sketch of the wellposedness result for (which is well 

known, see [6] and [3]) by introducing another approximation scheme. This scheme becomes 
also very useful to prove an I/°°-estimate for the approximation of dI^-(^ (independent of the 
approximation parameter) via the maximum pronciple. This kind of estimate at the e-level will 
guarantee that the same bound remains valid also for dl^yff. This fact will be rather crucial in 
proving (|1.5p . 

Under the assumption that the estimate (jl.Sp for the Lagrange multipliers is known, 
in Section [3] we give the conditional proof of the existence of an exponential attractor for the 
solution semigroup S{t) associated with equation (jl.ip . This result is obtained using some 
modification of the so-called method of Z-trajectories which was originally introduced by Malek 
and Necas in [H] and is widely used nowadays in the attractors theory, see [101 US] and references 
therein. 

Section[3]is the key part of the paper and it is devoted to verifying estimate (jl.Sp for the 
case in which the boundary dJ(f of the convex set J(f is regular enough. The proof is based on 
the maximum principle and the Kato inequality together with rather delicate construction of 
the approximating potentials. Unfortunately, the most relevant (from the applications point of 
view) choice of J(f is that of the simplex p.3p which does not fit with the regularity assumptions 
on the boundary of J(f . Thus, we have to consider this particular case separately using a different 
approximation scheme for the subdiffcrcntial dljf^. As a result, the existence of an exponential 
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attractor is proved for the case in which the boundary of is smooth and for the particular 
choice of Jif as a simplex (jl.3|) (although we expect that this result should be true for any closed 
and bounded convex set The section closes with some discussion on possible extensions of 
our result to more general classes of reaction-diffusion equations. 

Finally, in the last Section^ we study the convergence of the exponential attractors M.^ 
for the regular approximating problems to the limit exponential attractor oi the singular 
problem For simplicity, we restrict ourselves there only to the case of the simplex p.3p 

and verify that 

dist(r^(o))"(-^--^o) <C£«, 

where dist*'''" stands for the symmetric Hausdorff distance and the positive constants C and k 
are independent of e. The result is based on the proper application of the abstract theorem on 
perturbations of exponential attractors proved in [M] to the obstacle problem 

2 Well-posedness and regularity 

The aim of this section is to recall some known facts about the solutions of reaction-diffusion 
problem p.ip with obstacle potential and to formulate some additional estimates which will be 
crucial for what follows. We start with 

A word on the notation: The unknown function u is actually a vector valued function but, 
for the sake of simplicity, will be denoted as a scalar valued function. Consequently, also the 
functional spaces we will use in the course of the paper we will have a "scalar" notation. This 
means that a notation like, e.g., I?' will be preferred to a (more precise) notation like (L^(£l))". 
The same applies to dualities ((•, •)) and scalar products ((•, •)). Moreover, we will indicate with 
same symbols Jif and Jjf the convex in R" and its indicator function and their realization in 
!?• (see the next proposition 12. 4p . Thus, the definition of (weak) solutions of our problem is. 

Definition 2.1. A hmction u — u{t, x) is a solution of the obstacle problem if u{t, x) £ 
for almost all (t, x) G [0, T] x n, 

ueC([0,r];L2)nL2(0,T;i/i), dtu e L^{0,T; H'^), (2.1) 

and the following variational inequality holds for almost every t G (0, T] 

{dtu{t), u{t) - z) + (Vu, V(u - z)) < A(u, u - z), for any z e X. (2.2) 

The next theorem is a standard result in the theory of the evolution equations associated 
with maximal monotone operators (see the seminal references [S|, [B] and [3]). 

Theorem 2.2. [Well posedness] Let be a closed and bounded convex set containing the 
origin in R". The, for any given measurable uq such that uq{x) G .J(^ for almost all x € 
there exists a unique global solution u of problem (jl.ip in the sense of definition \2.1\ Moreover, 
u{i) G and dtu(t) G for t > and the following estimate holds: 

\\dMt)\\L- + \\uit)\\H2<C^^, (2.3) 

where the constant C depends on A and ,J€ , but is independent of uq. If, in addition, uq G 
n Hq , then the term with in the right-hand side of (j2.3p can be removed. 
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Proof. Although this result is well-known, for the convenience of the reader, we briefly 
recall one of its possible proof and deduce the regularity estimate (|2.3p . 

We first note that the uniqueness follows immediately from the variational inequality 
(|2.2p . Indeed, let u and v be two solutions of (jl.ip . Then, taking a sum of (|2.2p for u{t) and 
z = v{t) with ()2.2p for with z = we have 

{dtu{t) - dtv{t),u{t) - v{t)) + \\V,{u{i) - v{t))\\l, < \\\u{t) - v{t)\\l,. 

Denoting now w{t) :— u{t) — v{t) we arrive at the differential inequality 

1 d 



2 dt 



MMh + \\^Mmh < MHmh- (2.4) 



Applying the Gronwall inequality and using that, by definition u,v G C'^([0, T]; i^) for any 
T > 0, we finally have 

Wwimh + f \\Hs)fm ds < Ce'^\\wm\%, (2.5) 
Jo 

where the constant C is independent of u and v. Thus, the uniqueness holds. 

The existence part is slightly more delicate and requires the approximations of the 
singular convex potential Ij(f by suitable regular ones i^e . In fact, thanks to the uniqueness result 
this approximation can be done in several ways. The typical choice is usually the Morcau-Yosida 
approximation (see [Hj). However, in view of the next results (see in particular Proposition [^3] 
and Theorems 14.11 and 14. 6p , it is more convenient to adopt another kind of approximation. To 
this end, we let M{u) be the distance from the point u G M" to the convex set . namely the 
real valued function M 

Af(u) := dist(u,^). (2.6) 

Then, M is convex, globally Lipschitz continuous and smooth outside J(f . In addition there 
holds that 

M{u) > 0, M{u) = 0, if u e Jf, (2.7a) 

\VM{u)\ = 1, if u ^ J^. (2.7b) 

Now for any £ > 0, we introduce the real function 

- e)^, z > e. 

Finally, the desired approximation is defined a 

F,{u) feiM{u)). (2.9) 

Then, obviously, F^{u) is convex and smooth (at least C^'^) and F^{u) = F^{u) = for all 
u G J^, where F^ denotes the gradient of F^. Wc thus consider the following approximation of 
the problem p.ip : 

idtu- A^u + F^iu) = Xu, ^2^^^ 
where the initial data uq G L°°{fl) and uo{x) e JXf a.e. in O. 
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Our next step is to obtain a number of uniform (with respect to £ — *■ 0) estimates for 
the sohitions of the regular problems (|2.10p . We start with the usual L^-estimate. Indeed, 
multiplying equation (j2.10p by u, integrating over 17 and using the obvious facts that 

\{F',{u). - > -C, {F^u), u) = (F^iu) ~ F^O), u) > 

where the constant C depends only on A and on Q but is independent of e \ 0, we arrive at 

ljtMi)\\h + W^Mmh + liK(.^^),u) < c. 

Applying the Gronwall inequality to this relation, we see that 

\Wmh + j'^' \\u{sWhI + {KHs))Ms))ds < C. (2.11) 

Here we have implicitly used that .J(f' is bounded and, therefore, ||wo|ji2 is uniformly bounded 
for all admissible initial data uq. 

Next, we obtain the uniform energy estimate for the solutions of ()2.10p . To this end, 
we multiply the equation (|2.10p by dtu(t) and integrate over x: 

|(l|v.u(i)||i. + 2(^^,(u(t)),i)) + < A|hi(t)||i.. 

Integrating this estimate with respect to time, using (12.1111 together with the obvious estimate 

0<{F,{u),l)<{FUu),u) 

(which is an immediate consequence of the convexity of F^) and arguing in a standard way (see, 
e.g., [29]), we deduce the desired TJ^-energy estimate 

lkW!l?fi+(i^e("(i)),l)+y[*^Vt«(s)|li.ds<C^, t>0, (2.12) 

where the constant C is independent of uq and s. 

We are now ready to verify the analogue of (|2.3p for the approximate solutions. We 
start with the estimate of dtu{t). To this end, we differentiate (|2.10|) with respect to time and 
denote v{t) := dtu{t). Then, we have 

dtv - A^v + F^'{u)v ^ Xdtu. 

Multiplying this equation by v{t) integrating with respect to time, using (j2.12p together with 
the monotonicity of F^(u) and arguing in a standard way (see again |29j for details), we have 

WdMml. < C^^±P^ (2.13) 

for some constant C independent of uq and e. 

Now, in order to obtain the iJ'^-part of estimate (|2.3p , it is sufficient to rewrite equation 
(|2.10p in the form of elliptic equation (for every fixed t) 

AMt)-F^{u{t))^dMt)-Xuit), M(t)| =0, (2.14) 
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multiply it by Aa:u{t), integrate over fl and use the fact that, thanks to the convexity of F^, 
the term {F!,{u), —A^u) is nonncgativc. Then, due to (|2.13p . (|2.12p and the elliptic regularity 
result for the Laplacian (recall that the domain Q is assumed to be smooth enough), 

\\u{t)\\H2 < C\\AMt)\\L- < C{\\dtu{t)\\L2 + X\\u{t)\\L'2) < C'^, (2.15) 

where the constants C and C are independent of e and uq. Furthermore, from equation ()2.14p . 
we now conclude that 

Finally, we are now ready to pass to the limit e — > in equations (j2.10p and verify the existence 
of a solution u for the limit obstacle problem. Indeed, let e„ > be a sequence of positive 
numbers such that £„ — > as n ^ oo. Assume for the first that the initial data uq is smooth 
enough and verify the constraint, namely 

uo e i?^ n iJg and uoix) G ,Jf, a.e. in € fi. 

Let Un{t) = (i) be the corresponding solution of the approximating problem (j2.10p . Then, 
analogously to the derivation of (|2.3p for u„, but using in addition that uq is regular, we have 

WdtUnimL^ + WUnimH- + 1 1 Ju„ (<)) 1 1 + 1 1 (u„ (0) 1 1 < Q{\\uo\\h-), (2.16) 

where the function Q is independent of n and t > 0. Thus, up to not relabeled subsequence, 
we have that m„ converges weakly star in 

X L°° (0, T; H^) n W^-^{0, T; L^) 

and strongly (thanks to [13 Cor. 4]) in C°([0, T];H^) to some function u e X D C°([0, T];H^). 
So, we only need to check that u is a desired solution of p.ip in the sense of Definition 12.11 
Indeed, the regularity of u is obvious. The fact that u{t, x) G J(f for almost {t,x) £ [0,T] x ft 
follows in a standard way from the uniform bounds for F^^(un{t)) in the L^-norm and from 
the fact that F^^ (w) oo for all w ^ J^. So, we only need to check the variational inequality 
(|2.2p . Let z = z{x) be any admissible test function. Then, using the monotonicity of F'^{u) and 
the fact that F'^(v) = for all v £ , we see that 

{F',{u{t)),u{t) - ^) - {KHt)) - Ki^),^ -z)>Q. 

Multiplying now ()2.10p by u(t) — integrating over Q and using the last inequality, we arrive 
at 

(9tU„(t),U„(t) - z) + {V xUn{t) ,V x{Un[t) - z)) < A(u„ (t) , U„ (f) - z) . 

Passing to the limit n ^ oo in this inequality and using the above convergences, we see that 
(|2.2p holds for almost any t and all admissible test functions z. Thus, we have proved the 
existence of the desired solution u for the case of smooth initial data uq. 

For general initial data mq, the corresponding solution can be now constructed by 
approximating the non-smooth initial data uq by smooth ones Uq and passing to the limit 
n — > oo. To conclude, we finally recall that the continuity of the limit function near zero {u e 
C([0, r], L-^)) follows from the global Lipschitz continuity (j2.5p ). Theorem 12.21 is proved. I 

Remark 2.3. The above proof reveals that the concrete form of the functionals Fg{u) is not 
essential for the proof of the existence result. In fact, the solutions of the approximating 
problems (|2.10p will converge to the unique solution u of the obstacle problem if the following 
assumptions are satisfied: 
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1. The functions i^^ are convex and regular enough 

2. i^,' ^ as e \ 0. 

Where the last condition means that 

i) F^{v) 0, Vw e JT, 

In other words, we only need the subdifferential F^{u) of F^ to be an approximation, via graph 
convergence, of the subdifferential (see, e.g., [U Proposition 3.60 and Theorem 3.66]). 

However, as we have already mentioned, the choice of suitable approximating functionals -F'e(w) 
turns out to be crucial for our method of estimating the dimension of the global attractor. Thus 
we will construct this approximation in a rather specific way depending on the structure of the 
convex set (see, however, the next sections). 

For the long time analysis it is actually more convenient to reformulate the variational inequality 
p.2p in terms of an equation coupled with a differential inclusion for the subdifferential of the 
indicator function of (actually of its realization) . This subdifferential, which a priori should 
be undestood with respect to the duality Hq-H^^, actually makes sense in and thus almost 
everywhere in U (see Theorem 12.21 and the next Proposition 12. 4[) . Thus, we introduce the 
function (named Lagrange multiplier in what follows) 

hu{t) -.^ -dtu(t) + AMt) + Xu, Kit) e dlj(r{u). (2.17) 

Then, due to Theorem /i„ G L°°{t,T;L'^), for any r > 0. Moreover, the definition of 
subdifferential (w.r.t. the scalar product) gives (|2.2p . 

{hu{t),u{t) - z) > 0, for almost all t e [0, T] (2.18) 

and any admissible test function z = z{x). The last inequality can be also written in a point- wise 
form. 

Proposition 2.4. Let the assumptions of Theorem \2.2\ hold and let = hu(t,x) be the La- 
grange multiplier associated with the solution u{t) of problem (jl.ip . Then 

{hu{t,x),u{t,x) - Z)u'^ >0, (2.19) 

and almost all (t,x) G [0,T] x D,, which means 

hu{t,x) G dLx{u{t,x)), a.e. in [0,T] x Vt. 

Proof. Indeed, arguing as in the proof of Theorem 12.21 we see that 

/i„„(t) :=F,'Ju„(t)) hu{t) 

weakly-star in the space L°°{t,T; L^), for any r > 0. Thus, it is sufficient to verify f2.19p for 
the functions h^^ only. But these inequalities are immediate due to the monotonicity of F^ and 
the fact that F^{Z) = if Z G JT. Indeed, 

{hu^{t,x),Un{t,x) - Z)r^ = {F^^{Un{t,x)) - F^ ^{Z)) , Un{t , x) ~ Z)r^ > 0. 

Passing to the limit n ^ oo in that inequalities, we deduce (|2.19|1 and finish the proof of the 
proposition. I 
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Thus, the obstacle problem (jl.ip can be rewritten in terms of functions u and hu as follows: 
dtu — A^u + hu = Xu, in the sense of distributions, 

huit,x) e dlj(f{u{t,x)), for almost aU (i, x) e M+ X rj, (2.20) 

The next proposition shows that the function is, in a fact, globally bounded in the L°°-norm. 

Proposition 2.5. Let the assumptions of Theorem \2.2\ hold and let hu{t) be the Lagrange 
multiplier associated with the solution u{t) of problem p.ip . Then, h^ G L°°(M+ x Q) and 

\KmL^<C, (2.21) 

where the constant C depends only on and A (and is independent of u and t>0). 

Proof. Thanks to the lower semicontinuity of norms with respect to the weak convergence, 
we will verify ()2.21|) only in the case in which /i„ in equation (|2.20p is replaced by its approxi- 
mation hu„{t) :— F^^{un{t)). To this end, we test equation (|2.10p in the scalar product of M" 
with VM(m) (where the function M is the same as in the proof of Theorem l2.2p . and use that, 
due to the convexity of M, 

(-A^u, VM(u))r,. = -A^(Af(u)) + (7^(A/(u))V:,u, V^u)r,. > -A^(Af(u)) 

where TC{AL) denotes the Hessian matrix of M (actually, H{M){u) does not exist if u € dJXf, 
but the inequality still holds and can be easily verified, say, by approximating the non-smooth 
convex function M by the smooth convex ones). Then, we arrive at the differential inequality 
for the function M{u) 

dtM{u) - A^M(u) + fi{M(u)) < X{u,VAI{u)), M{u)l^^ = 0, 

where we have implicitly used that |VA/(u)p = 1 for u ^ J(^ (see (|2.7bp ). Furthermore, since 
|u| < M(u) + C,jfr, where C,jfr := diam(J^), we finally have 

dtMiu)-A,M{u) + f^jMiu))<XM{u) + XCK, M{u)l^^^O. (2.22) 

Applying the comparison principle to the scalar parabolic equation (j2.22p . we see that 

M{u{t,x)) < Ve„, for a.e. {t,x) e M+ x f}, 

where the constant (w.r.t. x and t) v^^ > is the solution of the following equation 

In addition, from (|2.8p . we see that ^ as £„ (i.e. n / -foo). For this reason 

/^jAf(u)) < /;jt;,J = Xv,„ + XCx < Cx. 

and, therefore, 

l^n„(t)| =/e„(M("«(t)) <Ca 

uniformly with respect to n / oo. Finally, passing to the limit n /* oo, we arrive at (|2.2ip 
(with C ~ Adiam(^)) and finish the proof of the proposition. I 
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Remark 2.6. In contrast to Theorem 12.21 and Proposition [23] which are based only on the en- 
ergy type estimates which are vahd for much more general equations, e.g., with non-scalar 
diffusion matrix, etc., the L°°-estimate obtained in Proposition 12.51 is based on the maxi- 
mum/comparison principle and requires the diffusion matrix to be scalar. In particular, we 
do not know whether or not this estimate remains true even for the case of diagonal, but 
non-scalar diffusion matrix. 

As direct consequence of the previous Proposition, we have that (|2.20p could be understood as 
the heat equation 

dfU — AxU — Xu — hu 

with the external forces belonging to . Thus, the parabolic interior regularity estimates give 
(see. e.g., [H]) 

Corollary 2.7. Let the assumptions of Theorem ()2.2p hold and let u{t) be a solution of problem 
(|l.ip . Then, for every ly > 0, u{t) £ C'^~'^(fl) for t > and the following estimate holds: 

Mt)\\c^-^<C.^^, (2.23) 

where the positive constants and a are independent of t and u. If, in addition, Uq € 
C^-" n i?o\ then ^^n^ holds with a = 0. 



3 Global and exponential attractors 

The aim of this section is to study the long-time behavior of solutions of problem p.ip in terms 
of global and exponential attractors. We first recall that, due to Theorem 12.21 equation p.ip 
generate a (dissipative) semigroup {S{t), t > 0} in the phase space 

$ = := {m e : u{x) G for almost all x G n}, (3.1) 

i.e., 

S{t):^^^, S{t)uo:^u{t), (3.2) 

where u{t) is the solution to (|l.ip at time t. Moreover, due to estimate ()2.5p this semigroup is 
globally Lipschitz continuous in the i^-metric 

\\u,it) ^ U2{t)\\l. + f \\u,{s) - U2{s)\\]j, ds < Ce^*||«i - ulWl,, (3.3) 
Jo 

where the positive constants C and ^ are independent of t and of the initial data Uq,Uq G $. 
In addition, due to Corollarv l2.71 we have the following regularization estimate 

WdMmL^ + hWllc-^(O) < G.^^, (3.4) 

where v > is arbitray and the positive constants and M are independent of the initial 
condition uq and of t. These two estimates immediately imply the existence of a global attrac- 
tor A for the semigroup S{t) associated with the obstacle problem (jl.ip . We recall that, by 
definition, a set ^ C $ is a global attractor for the semigroup S{t) : <& ^ $ if 

1. The set A is compact in $; 
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2. It is strictly invariant: S{t)A ^ A, t > 0; 

3. For every neighborhood O — 0{A) of A in $ there exists a time T — T{0) such that 

S{t)<S> C ©(y^), t > T. (3.5) 

Remark 3.1. The attraction property (j3.5p is usually formulated not for the whole phase 
space but for the hounded subsets of it only. In our case, however, the whole phase space 
$ is automatically bounded (since is bounded in M"), so we need not to use bounded sets 
to define the attractor. Note also that the attraction property p.5p can be reformulated as 
follows: 

distioc {S{t)<i>, A)^Q as t ^ oo, (3.6) 
where distv{X^Y) := sup inf dv{x,y) is the non-symmetric Hausdorff distance between the 
subsets X and Y of the metric space V. 

Theorem 3.2. [Global Attractor] Under the assumptions of Theorem \2.2l the semigroup S(t) 
associated with the obstacle equation (|l.ip possesses the global attractor A in ^ which is bounded 
in C^~'^{fl) for every u > 0. This attractor is generated by all the trajectories of the semigroup 
S{t) defined for all t G R: 

A = Kl^^, (3.7) 
where K C L°°(K., <i>) is a set of all solutions of (|1.1|) defined for all t G M. 

Proof. The proof of this theorem is standard and well known (see, for instance [2], [H], 
or [H]). Indeed, due to estimate p.4p the semigroup S{t) possesses an absorbing set which 
is bounded in C^~'^(r2) and, therefore, compact in L°°(r2) and estimate (|3.3|) guarantees that 
the semigroup has a closed graph. Thus, all assertions of the theorem follow from the abstract 
attractor existence criterium, see [5]. I 

Recall that the semigroup S{t) associated with the obstacle problem ()l.ip possesses a global 
Lyapunov function of the form 

C{u) -.= 11^ :M\h-X\\u\\l,. (3.8) 

Indeed, using the test function z = u{t — h) in the variational inequality (|2.2p . dividing it by 
h > Q and passing to the limit ^ 0, we arrive at 

\\dtu{t)\\l. + j^C{u{t))<Q. 

Therefore, according to the general theory (see e.g., [2]), every trajectory u{t) = S{t)uo tends 
as t ^ oo to the set TZ of all equilibria of problem (jl.ip 

dist{S{t)uo,TZ) ^ 0, as t ^ oo. 

However, in contrast to the case of regular systems, the equilibria set TZ is generically not discrete 
for the obstacle type singular problems. Thus, in our situation, the existence of a Lyapunov 
function does not allow to obtain the stabilization of every trajectory to a single equilibrium 
even in "generic" situation. In addition, the semigroup S(t) is not differentiable with respect to 
the initial data (it is in fact only globally Lipschitz continuous), so we are not able to construct 
the stable/unstable manifolds associated with an equilibrium. Thus, the so-called theory of 
regular attractors is not applicable to equations of the type p.l|) . Moreover, due to the above 
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mentioned non-diffcrcntiability, the standard way of proving the finite-dimensionality of the 
global attractor based on the volume contraction method does not work here. So, the existence 
of the finite-dimensional reduction for the associated long time dynamics becomes a non-trivial 
problem which, to the best of our knowledge, it has not been yet tackled. In this paper we will 
prove that the global attractor for (jl.ip has finite fractal dimension by using the concept of the 
so-called exponential attractor and the estimation of the dimension based on the proper chosen 
squeezing/smoothing property for the difference between two solutions. This method has the 
advantage that it does not require the differentiability with respect to the initial data The 
existence of an exponential attractor is interest in itself. In fact, we recall once more that the 
global attractor represents the first (although extremely important) step in the understanding 
of the long-time dynamics of a given evolutive process. However, it may also present some severe 
drawbacks. Indeed, as simple examples show, the rate of convergence to the global attractor 
may be arbitrarily slow. This fact makes the global attractor very sensitive to perturbations 
and to numerical approximation. In addition, it is usually extremely difficult to estimate the 
rate of convergence to the global attractor and to express it in terms of the physical parameters 
of the system. In particular, it may even be reduced to a single point, thus failing in capturing 
the very rich and most interesting transient behavior of the system considered. The simplest 
example of such a system is the following ID real Ginzburg-Landau equation 

dtu = ed^u + u — u^, x G [0, 1], = = —1. 

In that case, the global attractor A = { — 1} is trivial for all e > 0. However, this attractor is, 
factually, invisible and unreachable if e is small enough since the transient structures (which 
are very far from the attractor) have an extremely large lifetime T ~ e^/ ^ . 

In order to overcome these drawbacks, the concept of exponential attractor has then 
been proposed in [5]) to possibly overcome this difficulty. We recall below the definition of an 
exponential attractor adopted for our case, see e.g., [8] and [23] for more detailed exposition. 

Definition 3.3. A compact subset Ai of the phase space $ is called an exponential attractor 
for the semigroup S{t) if the following conditions are satisfied: 

(El) The set M is positively invariant, i.e., S{t)M C M for all t > 0; 
(E2) The fractal dimension (see, e.g., [HIIIH]) of in <I> is finite; 

(E3) The set M attracts exponentially fast the image the phase space Namely, there exist 
C,P > such that 

distL^{S{t)^,M) <Ce-'^\ Vt > 0. (3.9) 

Thanks to the control of the convergence rate {E3) it follows that, compared to the 
global attractor, an exponential attractor is much more robust to perturbation (usually it is 
Holder continuous with respect to the perturbation parameter, see Section [5] below) . However, 
since the the exponential attractor M is only positively invariant (see (El)), it is obviously 
not unique. Thus, the concrete choice of an exponential attractor and its explicit construction 
becomes essential. We recall also that, in the original paper [8] the construction was extremely 
implicit (involving the Zorn lemma) and this fact did not allow to develop a reasonable pertur- 
bation theory. This drawback has been overcome later in [5] and |10j where an alternative and 
relatively simple and explicit construction for the exponentially attractor has been suggested. 
Note also that the construction of [TU] gives an exponential attractor which is automatically 
Holder continuous with respect to the reasonable perturbations of the semigroup considered and 
this somehow resolves the non-uniqueness problem. We refer the reader to the recent survey 
P5] for the detailed informations on the exponential attractors theory. 
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The proof of the existence of an exponential attractor for our obstacle problem (jl.ip 
will be organized as follows. First of all, we establish (in a quite standard way) the existence 
of an exponential attractor under a crucial additional (for this moment only) assumption on 
the differences of the Lagrange multiplies h^^ and of two different solutions ui{t) and U2{t) 
of (jl.ip . This assumption is the core of the exponential attractor existence Theorem and it is 
the main result of the paper. The next section [4] will be dedicated to its proof in two different 
situations: the convex set J(f is smooth, or the convex set is an n-dimensional simplex. 

Assumption There exist a closed positively invariant absorbing set Bq d for the 

semigroup S{t) such that, for every two solutions ui and U2 of ()2.20p starting from Bq (i.e., 
Ui{0) e Bq), the corresponding Lagrange multipliers /i^ and h^^ satisfy the following estimates: 

\\hu, - /lnJ|LM[0,l]xO) < •^llui(O) - U2(0)||l2, (3.10) 

where the constant is independent of u\ and U2 . 

Under Assumption in the next Theorem we prove the existence of an exponential 
attractor (and thus the finite dimensionality of the global attractor) . 

Theorem 3.4. Let the assumption of Theorem \2.^ hold and let, in addition, the semigroup 
S(t) associated with equation (jl.ip satisfy Assumption .f£ . Then, S{t) possesses an exponential 
attractor A4 CZ ^ in the sense of Definition \3.3[ Moreover, A4 is a bounded subset of C^~^ , 
for all V > Q. Finally, the global attractor A constructed in the Theorem \S.2\ has finite fractal 
dimension. 

Proof. Recall that, since Bq is a semi-invariant absorbing set for S{t), it is sufficient to 
verify the existence of the exponential attractor for the restriction of S{t) on Bq only. As 
usual, we first verify the existence of such attractor for the discrete semigroup generated by the 
map S = S{1) and then extend to the continuous time. To this end, we will use the following 
abstract exponential attractor existence theorem suggested in 0. 

Lemma 3.5. Let Jif and J^i be two Banach spaces such that is compactly embedded to H . 
Let Bq be a bounded closed subset of Jf° and a map 

S : IBo — >■ IBq 

be such that 

- §&2||^i < K\\bi - b2\\.^, 6i, 62 e Bo, (3.11) 

where the constant K is independent of bi and 62. Then, the discrete semigroup {§(«), n G N} 
generated on Bq by the iterations of the map § possesses an exponential attractor, i.e., there 
exists a compact set Add C Bq such that 

(El) Aid is positively invariant: EAAd C A4d; 

(El) The fractal dimension of Aid in is finite: 

dim/(7Wd,.^) < A/ < 00; 

and 

(E3) Aid attracts exponentially the images o/Bq under the iterations of the map Mq: 

dist^(§(7i)Bo,Xd) < Ce-''". 

Moreover, the positive constants M , C and k can be expressed explicitly in terms of the squeezing 
constant K , the size of the set Bq and the entropy of the compact embedding C . 
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We will use the so-called method of "^-trajectories (introduced by Malek and Necas in 
[12], sec also and [25]) in order to construct the proper spaces and J^i and to verify the 
assumptions of Lemma 13.51 

Namely, let us consider the trajectory space Bq consisting of the pieces of trajectories 
of the solution semigroup S(t) of length one starting from Bq: 

Bo := {u e L°°([0, 1], $), u(0) = uq e Bo, u{t) = S{t)uo, t £ [0, 1]}. (3.12) 

Then, there is a one-to-one correspondence between i?o and Bo generated by the solution map 

T : Bo^Bo, (Tuo)(i) S{t)uo 

and, therefore, we may lift the semigroup S{t) : Bq ^ Bq to the conjugated semigroup §(t) 
acting on the trajectory space Bo: 

§(i):Bo^Bo, §(t) :=To5(t)oT-^ (3.13) 

We intend to apply Lemma [33] to the map § — §(1) acting on the trajectory space Bo. To this 
end, we define the spaces and as follows: 

Jf' := L2(0,1;L2), J^i := ^^(0, 1; iJg^) n W^^'1(0, 1; i/""), (3.14) 

where s > max{l, N/2} is a fixed exponent. Then, obviously, the embedding C Jif is com- 
pact and we only need to check the smoothing property (|3.1ip . To this end, we need (together 
with Assumption ^ and estimate p.3p ) the following additional regularization property of the 
semigroup S{t). 

Lemma 3.6. Letui{t) andu2{t) be two solutions of problem (|l.ip . Then, the following estimate 
holds: ^ 

\\ui{l)-U2{l)\\l2<{2X + l) [ \\ui{t) ^ U2{t)\\l2 dt. (3.15) 

Jo 

Indeed, multiplying the differential inequality ()2.4p by t and integrating t G [0,1], we 
arrive at p.l5p . 

We are now ready to verify the smoothing property (|3.1ip . 

Lemma 3.7. Let the assumptions of Theorem \2.Si hold and, in addition, Assumption ^ be 
satisfied. Then the following estimate holds for every two solutions ui{t) and U2(t) such that 
Wq e Bq, i ^ 1,2; 

- W2|1l2(1,2;H1) + \\dtUl ~ dtU2\\L^l,2;H--) < MWl - "2 |1 L2 (oa;L2) , (3.16) 

where the constant L is independent of ui and U2. 
Proof. First of aU we recall that, for v G L^(l, 2; iJ^"*), 



m\mi,2;H-) = sup 



2 



(v, ip)dr 



where the sup is taken over the ip G L°°{1,2; Hq) such that ||</'||L=°(i,2;_f/g) — 1 a-ud the duality 
pairing is of course between H~'' and Hq. Consequently, there holds 

2 

\\dtui{t) - dtU2{t)\\H-sdt 

<l \\ui{t)-U2{t)\\Hidt+ f \\hu,{t) - hu,{t)\\L^dt + \ f \\ui{t)-U2{t)\\L2dt 
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(here we have impHcitly used that s > max{l, iV/2} which impUcs that Hq C L°°). Using now 
Assumption together with the global Lipschitz continuity (|3.3p . we have 

- "2i|L2(l,2;ffl) + \\9tUl - 9fli2||Ll(l,2;ff-) < C'll'"l(l) ^ ■"2(1)|1l2- 

Combining this estimate with p.lSp . we arrive at (j3.16p and finish the proof of the lemma. I 

We are now ready to finish the proof of Theorem 13.41 Indeed, we have verified that the map 
S := §(1) satisfies all of the assumptions of the abstract Lemma l531 and. therefore, the semigroup 
{§(n), n € N} possesses an exponential attractor in the trajectory space Bq endowed with 
the topology of =^ = i^(0, 1; L^). Projecting it back to the phase space Bq C ^ via 

Md:^Md\^^^ (3.17) 

and using estimate (|3.15p . we see that Md is indeed the exponential attractor for the discrete 
semigroup {S{n), n £ I,} acting on Bq (endowed with the topology of L^). In addition, we see 

that 

MdC S{l)Bo C Sil)<P (3.18) 

which implies that ||9tu(i)||/^2 < C for every trajectory u{t) starting from Aid (due to estimate 
(|2.3|) ). Thus, thanks to p.3p . the map (t, uo) — * S{t)uo is globally Lipschitz continuous on 
[0,1] X Md (in the R x L^(ri)-metric). Thus, the desired exponential attractor for the solu- 
tion semigroup S{t) with continuous time can be now constructed by the following standard 
expression (see [5] for the details): 

M := Ute[o.i]Md- 

Note that; up to now, we have verified that Ai has the finite fractal dimension and possesses 
the exponential attraction property in the topology of only. In order to verify that these 
properties actually hold in the _L°°-topology of the phase space we use (in a standard way) the 
additional regularity of A4 and a proper interpolation inequality. Indeed, by our construction, 
A4 C 5'(1)$ and, therefore, due to estimate (|2.23p . A4 is globally bounded in C^"". Using now 
the following interpolation inequality 

\\u-v\\l^ <C,\\u~vrL^\\u~v\\l7,^^, 0<K<1, (3.19) 

we see that the dimension of Ai is finite not only in , but also in and that the attraction 
property holds in L°° as well. Thus, the desired exponential attractor M in the phase space $ 
is constructed and Theorem 13.41 is proved. I 

Remark 3.8. Note once more that the method of proving the conditional result of Theorem 
13.41 (which is more or less standard variation of the Z-trajectories method) is widely used now- 
a-days in the exponential attractors theory, see e.g., [13] and the references therein. Thus, the 
major difficulty here (and the major novelty of the paper) is related with the verification of the 
Assumption J^f for the solution semigroup S{t) associated with the obstacle problem (|1.1|1 . 

4 Estimates on the difference of the Lagrange multiphers 
and verification of Assumption =Sf 

The aim of this section is to verify the Assumption ^ under some additional assumptions on 
the structure of the convex set We start with the case in which ^ has a smooth boundary. 
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4.1 The case of regular Jif . 

The main result of this subsection is the foUowing Theorem. 

Theorem 4.1. Let the assumptions of Theorem \ 2.2\ hold and let, in addition, the boundary 
S = he smooth enough (at least, C^'^j. Then, the solution semigroup S{t) associated with 
equation (jl.ip satisfies Assumption ^ . 

Proof. The proof is based on an argument that combines a proper choice of an approxima- 
tion scheme and the maximum principle similar to the one devised to prove Proposition 12.51 
However, since the function M(u) = dist(u, J?r) used to define in the existence Theorem l2.2l 
is not smooth enough near the boundary, we should introduce another approximation of the 
singular potential. The smooth correction of M is given by the following Lemma. 

Lemma 4.2. Let .J^f he a hounded convex set with the C^'^ -smooth boundary S. Then, there 
exists a function M : K" M with the following properties: 

MeC^'^W)- (4.1a) 

M is convex; (4- lb) 

S = d.je = {z e M" such that M{z) = 0} ; (4.1c) 

\VM{z)\=e{M{z)), (4.1d) 

where ■& = d{z) is a monotone increasing function which is smooth near z = and such that 
0(0) ^ 0. 

Proof. The function M can be constructed as follows. First of all, for any 5 > 0, we 
introduce the following set: 



S-f, := 



[z^ek: dist(zo,^) =5}. (4.2) 



Then, we consider the subset J(f-s C ^ whose boundary is S^s- It turns out that, being 
convex, the domain J(f-j is convex too. Moreover, by possibly taking 5 small, the boundary 
S-s has the same regularity of S. The candidate for M is thus the function 

M{z) := [dist(z, ,je-s)f - 5^. (4.3) 

In fact, M is clearly convex and regular. Moreover, by choosing d = ■d[w) = 3(w + and 
using the identity 

|v(dist(z,jr_5))| = 1, 

one sees that the condition (|4.1dp is satisfied. 

Finally, by possibly taking a small 5, the smoothness of the boundary S entails the 
validity of (|4.1cp . Note also that M has the additional property of qualifying the fact that 
V £ Jf, namely it turns out that v E ^ if and only if —8^ < M{v) < 0. I 

Now, we introduce the approximations [u) of the indicator function Ijff as in (j2.9p using the 
above defined function M. Let now u{t) and v{t) be two solutions of the singular equation 
(|2.20p with initial conditions 



' 



with sufficiently small positive v (actually, = 1 is sufficient for what follows). Then, let Ue{t) 
and Ue(t) be their approximations, namely the solutions of the approximate problems (j2.10p 
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with initial data uq and vq. As one can easily check, Proposition 12 . 51 and Corollarv l2. 71 still hold 
for this new approximation entailing the validity of the following two estimates 

heWllc^- + \\Ve{t)\\c2-.^ < Qdlliollc^-^ + \\vo\\c^-.'), t > 0, (4.4) 

and 

\\F^{uS))U^+\\FUvS))U^<C, (4.5) 

where the monotone function Q and the constant C are both independent of s. Thus, up to 
not relabeled subsequence, 

weakly strar in L°°([0,T] x fi). Thus, in order to prove the theorem, we only need to verify 
that 

\\F'^{u^) - i^;(we)||Li([o,i]xn) < Q(ho||c^-- + ||wo||c2-OII"o - ^olU^ (4.6) 

with the function Q not depending on e and then by semicontinuity we obtain the desired 
estimate p.lOp in the limit e \ 0. This estimate will hold for every two trajectories u and v 
such that Uq and are bounded in C^~^ . In order to construct the desired positively invariant 
absorbing set Bq, it is sufficient to take a ball 

B := {u e C2-''/2 p ^1 p ||m||c.2-. < i?} 

with a sufRciently large radius i?. Then, thanks to Corollary 12.71 B will be an absorbing set 
for the semigroup S{t) associated with equation p.ip and, therefore, the closure Bi := of 
the set B in the topology of $ will be a bounded in C^^" and closed in $ absorbing set for the 
semigroup S{t). Applying CoroUarv 12 . 71 again, we see that the set 

Bo ■■= U Bi (4.7) 
t>o 

will be the desired positively invariant, closed in $ and bounded in C^~''{n) absorbing set for 
semigroup S{t) and estimate ()4.6p will guarantee that (j3.10p will hold uniformly with respect 
to all trajectories starting from Bq. 

Thus, it only remains to verify the uniform estimate (|4.6p . The first step is, roughly 
speaking, to reduce p.ip to an equation with scalar constraint (actually to an approximation 
of). This is done using the function M as we did in ()2.22p . Indeed, testing equation (j2.10p in 
the scalar product (•, •)!{" of R" with VA/(it) and using the condition (j4.1dp . we obtain 

dt{M{ue)) - A(M(Me)) + T',{M{u,)) + D{u,) = 0, (4.8) 

withP^{z) := f'e{z)9^{z) and 

N n 

^(") M'lk{u)d^^u,d^^Uk ~\Y,M'^{u)uj, (4.9) 

2—1 j,k—l j 

where {Mj^k}^ ^.^i denotes the entries of the Hessian matrix 'H{M) of the function M (and the 
analogous equation holds for u^). Let now := M{u^) — M{ve). Then 

dtw, - Aw, + [T'^iMiu,)) - T'^iMiv,))] + [D{u,) - D{v,)] = 0. (4.10) 
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We now multiply this equation by sgn and integrate over f2 obtaining 

= - [ (Diu,) - D{v,))sgnw,dx =: TZ, (4.11) 

where we have used the monotonicity of .^-"^(z) (recall in particular that i?(z) = {z + S^)"^^^) and 
the Kato inequality. To estimate the term with TZ we start to rewrite it in the more explicit 
form 

N 



N 
N 

+ / {'H{M{Ve))[dx,Ue - dx^Vs),dxiVs)^,,SgTL{We)dx 

+X {VM{u,),Ue)j^„sgii{w,)dx~ X {VM{v,),v,)g^„agn{we)dx. (4.12) 



Thus, thanks to the assumed C^'^ regularity on the function M (which follows from the smoth- 
ness of the boudary of Jf) and the C^-regularity of and which follows from (j4.4|) . it is 
not difficult to realize that 

mt)\\mn) < C\\u,{t) - veimm for any t > 0, (4.13) 

where the constant C depends on the -norms of uq and vq, but is independent of e > 

Integrating the differential inequality (|4.1ip with respect to t G [0, 1] and using ()4.13p 
and the analog of (jS.Sp for the solutions and of the approximate problems (|2.10p . we 
arrive at (recall also that 1X5(0) = uq and Ve(0) = vq) 

\\T^{M{ue)) - TUM{vMl^ dt < \\uo - voWl^ + C n{t)dt 



< \\u,iQ) - v,{0)\\li + C / \\u,it) - v,it)\\Hi dt 

Jo 

<Ci\\uo-vo\\l2, (4.14) 

where the constants C and Ci depend on the C^-norms of u{Q) and v{0), but arc independent 
of e. 

The next step is to estimate the difference f^{M{us.)) — f'^{M{v^)) in terms of the 
difference of T' {M{u^)) and JFg(Af (w^)). To this end, we note that fe{z) = if z < 0, so, 
without loss of generality, we may assume that, say, M{u^) > 0. Then, if -M(wj) < 

IfAMiu,)) ~ f,{M{v,))\ = fMiu,)) = e-^iMiu,))T',{Miu,)) 

< e-\Q)T',{M{u,)) ^ C\T',{M{u,)) - T',{v,)\. (4.15) 
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Let now M{v^) > 0. Then, 

\fUM{u,)) - .niMiv,))\ = \:FUM{u,))9-'{M{u,)) - ^;(A/K))0-2(Af K))| 

< e-'{M{u,)) ■ \J^UM{u,))~J'UM{ve))\+^'AMiv,))\e-'iMiu,)) - e-^{M{v,))\ 
< 9-'{0)\TUMiu,)) - TUM{v,))\ + \\F^{u,)\\l^9-H0)\9'{M{u,)) - ^^(MK))! 

< c(\TUM{u,)) - TUM{v,))\ + \u, ~v,\), (4.16) 

where the constant C is independent of s (here we have impHcitly used that the L°°-norms of 
F^{ue) are uniformly bounded, thanks to Thus, due to (|i7n|) . ((^7T5)) and (|i?T5)) 

' ||/^(Af(«e(t))) - /;(Af (z;aO))llLi dt < C\\uM ~ «,(0)I1l- 

Finally, since M is at least C^, 

||^^;(u,) - F;(iO||li = ||/^(Af(^,)A/'(^,) - /;(A/(i;,))A/'(iO||li 

< \\M'{u,)\\L^\\m-Hu,)) - /;(A//(z;,)|U. + ||/;(Af(^,))|U~||A/'(«,) - M'{v,)\\l. 

< CXmM{u,)) - f^iMivML^ + he - Veh^) (4.17) 

and 



\\F^{Miu,{t)) - F'{M{v,m\L. dt < C\\ue{0) - vMWl-- 
Thus, estimate (|4.6p is verified and Theorem 14. H is proved. 



Corollary 4.3. Let the assumptions of Theorem \2.2\ hold and let boundary S = dJ^ he of 
class C^'^ . Then, the semigroup S(t) associated with the obstacle problem p.ip possesses an 
exponential attractor M. in the sense of Definition in the phase space Moreover, the 
global attractor has finite fractal dimension. 

Indeed, this assertion is an immediate corollary of Theorems 13.41 and 14.11 



4.2 The case of an irregular convex J(f: the simplex (11.31) . 

As we saw in the former section, the the smoothness of the boundary S = dJ(f seems crucial 
for the method of verifying Assumption ^ suggested in the proof of Theorem l4.1l First of all, 
the C^'^ smoothness of M is necessary in order to obtain estimate (j4.13p . Secondly, the fact 
that 0(0) in (|4.1d[) . which is crucial to obtain estimates (|4.16p and (j4.15p ') implies that 
VA/"(z) ^ for all z € S. Thus, by the implicit function theorem, the boundary S also must 
be at least C^'^-smooth. However, from the possible applications to phase transition problems, 
one of the most important examples for the set J(f is the n-dimensional simplex, namely the 
set 

|p= (pi,...,p„) eR" such that < 1, > 0, i = 1, • • • , n| (4.18) 

which is a polyhedron and its boundary is only piece-wise smooth. Thus, Theorem 14.11 is 
not directly applicable here. Nevertheless, as we will show below. Assumption ^ remains 
valid for the non-regular case (|4.18p . Again, the verification of Assumption ^ will rely on an 
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approximation argument and, to this purpose, we will consider a slightly different approximation 
for the indicator function namely, let us introduce 

n n 

F,{u) ■.= MY^U,-l)+J2fe{-U^), (4.19) 

i=l i=l 

where the function is defined by (|2.8p . Obviously, F^{u) is convex and satisfies the conditions 
of Remark 12.31 Therefore, we can indeed use it for approximating the singular problem p.ip . 
Consequently, all of the estimates and convergences obtained in the proof of Theorem 1 2 . 21 hold 
for this new approximation. 

On the other hand, an inspection in the proof of the estimate in Proposition 12.51 
reveals that one of keys point was the structure of the approximations. It is also evident that 
the new defined approximation for the simplex is structurally different from the approximation 
introduced before; thus the L°° bound of Proposition 12.51 need a different proof. In particular, 
having in mind also the study of the approximation of the exponential attractor (see the next 
section), we formulate the analogue of Proposition 12.51 in a slightly stronger form by indicat- 
ing also the invariant regions at the e-lcvel. Thus, the global bound result on the Lagrange 
multipliers takes this form 

Theorem 4.4. Let e > be small enough and let the function Fg{u) be defined by (j4.19p . 
Then, there exists a positive constant p ( independent of e) such that the set 

je'e:={ue M", F^{u) < pe} (4.20) 

is an invariant region for the solution semigroup S'e(t) of the approximate problems (j2.10p . 
namely 

5e(t):$e^$e, t > 0, (4.21) 

where 

$e {w e L°°(17), u{x) £ X, for almost all X £ n}. (4.22) 

Finally, the approximations hu^{t) := F'^{u^{t)) (where u^{t) :~ S^{t)uQ, uq £ are uniformly 
bounded in the L°°-norm: 

\\huAt)\\L^m<C (4.23) 
where the constant C is independent of e and of the concrete choice of uq £ 

Proof. First of all, we need the following Lemma which clarifies the relations between 
|i^^(u)| and Fe(w). 

Lemma 4.5. Let the function F^ be defined by (|4.19p and the function fe{z) be given by ()2.8p . 
Then, there exist two positive constants ki and K2 ( independent of e) such that 

^F,{u) < \Fl{u)\^ < ^F,{u) (4.24) 

for all e > and all u £ . 

Proof. Indeed, let // -f',{-u,) and := f^{J2u, - 1). Then, 

n 71 n 

i—1 i—1 
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Let us consider two cases: 

Case I: /o = /e(E - 1) = 0- Then, (j^?^ simply reads as 

n 71 
i=l i=0 

Case II: /g = f'^(YliUi — 1) 7^ 0. Then, keeping in mind the definition of the function /e, 
we conclude that at least one of u^, i = 1, - • • ,n (say, m„ for definiteness) must be positive 
and, therefore, without loss of generality, we may assume that u„ > and, consequently, 
— — /^(— u„) = 0. Then, using the elementary inequality ab > —\{^ + ab"^), a > 0, we 
arrive at 

n~l n— 1 n— 1 

Z— 1 2 — 1 2—1 

Choosing the positive a = a{n) in an optimal way as a solution of the equation 

1 — a^^ ~ n{l — a) + a), 

we see that, in the second case 

n 

|F,'HP>0(n)El/;P 

2=0 

with 6{n) := ''+^-^/^^+^)^z ± > 0. Thus, in both cases 

n n 

|F;(u)P > 0(n) E = ^W(l/oP + E l/^l') • 

2=0 2=1 

Since the upper bound is obvious, we arrive at the following inequality 

Cn n \ n n 

E i/;(-"2)p + i/;(E"« - 1)1' < i^e'HP < E i/;(-^2)p+i/;(E"^'-i)i'- (4-26) 
2=1 2=1 / 2=1 2=1 

It only remains to note that, due to our choice ()2.8p of the function f^, 

\fUz)\' = ^Uz) 

and, consequently, (|4.26p implies (|4.24p and finishes the proof of the Lemma. I 

With the help of Lemma l4?5l it is now not difficult to complete the proof of Theorem l4.4l Indeed, 
testing equation (j2.10|l in the scalar product of M" with F^{u) and using that F'J{u) > 0, we 
have (compare with (|2.22|) ') 

dtF,{u) ~ A,{F,{u)) + \Fi{u)\^ < A(u,F;(u))k„. 

Using now Lemma 14.51 together with the obvious fact that 

\u\<F,{u) + C, 
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where C is independent of e ^ 0, we arrive at the differential inequaUty for the scalar function 
V{t) := F,{u{t)y. 

dtV-A,V+{^~X')V<Ci, V\g^^ = 0, (4.27) 

where the positive constants Ci and k are independent of e and the concrete choice of the 
solution u. Applying the comparison principle for the heat equations to (|4.27p . we see that the 
region {u, V{u) < p^} will be invariant if 

^ Ci 2Ci£ 

Pe > 



Thus, taking p := we see that the region ^ defined by (|4.20p will be indeed invariant with 
respect to the semi-flow Se{t) if e > is small enough. Thus, we only need to check estimate 
(|4.23p . To this end, it remains to observe that, due to Lemma 

e e 

Theorem 14.41 is proved. I 
We are now ready to formulate the main result of this subsection. 

Theorem 4.6. Let the assumptions of Theorem \4.4\ hold and e > be small enough. Then, 
for any two solutions Ue{t) :~ S^(t)uo and Vg{t) := Se(t)vo of the approximate problems (|2.10p 
starting from (uq,vq G ^e), the associated Lagrange multipliers hu^{t) := F^{u^{t)) and 
hv^it) := F^{ve{t)) satisfy the following estimate: 

1 

\\huAt)-h,At)\\LHn)dt<L\\uo-v4LHn), (4.28) 

where the constant L is independent of e, uq and vq. In particular, the limit senigroup S(t) 
associated with the obstacle problem (jl.ip satisfies Assumption I£ with Bq := <I>. 

Proof. In order to simplify the notations wc forget for a while the e dependence and will 
write u(t) := S^(t)uo and v{t) :~ S'j(<)ijo. Then, the function w{t) := nit) — v{t) solves 

n n 

dtw,- A,w, + [f^{-v,)^ f'^{~u,)] + [f'^CEu,-l)~ f^iY.V:-'^)] = >^w., z = ,n. (4.29) 

1=1 i=l 

Multiply now the ith equation of (|4.29p by sgnwi and integrate over Q. Then, taking the sum 
over i and using the Kato inequality, we arrive at 

n n n n 

9,iiHu^-Aii^iu.+^ii/;(-u.)-/;(-t'.)iiL^ <-([/;(E"^-i)-/^(E«^-i)]'E^g^^O- 

2—1 i—1 i—1 i—1 

(4.30) 

Let us estimate the right-hand side of this inequality. To this end, we first note that 

n 

I E I - 

1=1 
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In addition, in the worst case where the equahty holds, all Wi are of the same sign. Therefore, 
at that point the term ~ 1) ~ feiT^'i=i ~ 1) ^^^^ ^^^^ have the same sign and its 

product with sgnwi will be non-negative. Thus, 

n n n n n 

^[fUY.^^ - 1) - /;(E"' - 1)] -Jl'snw, < {n imj2^. - 1) - n{J2v^ - i)i 

z— 1 i=l i=l i— 1 i=l 

and 

n n n 

5t|klUi+E ll/e(-"0-./e(-t'.)llLi < \\\w\W^{n-\)\\r,i^U,-\)-re^V.,-\)\W. (4.31) 
i— 1 2—1 i—\ 

Next, in order to estimate the right-hand side of (|4.31|) . we sum all the equations (|4.29p and 
multiply the obtained relation by sgn(^"^j^ '"j^X]r=i ^O- Then, using again the Kato inequality, 
we have 

dAY^w.Aw + ^ii/;(E^. - 1) - /;(E«« - 

z— 1 i—\ z— 1 

n n 

< All E ^^\W + E - /^(-^OIIli. (4.32) 

1=1 i=l 

Multiplying inequality (j4.32p by and summing it to inequality ()4.3ip . we finally arrive at 

- / n n Ti \ 

+^ ii/^(E - 1) - /e(E - + E - rA-v^)\w 

\ 1=1 1=1 i=l / 

<A(^|HlL^ + ;^ilE'^^''ll^^)- (4.33) 
Applying the Gronwall inequality to this relation, we infer 

„1 n n n 

/ ii/e(E "»(^) - 1) - /e(E - i)iu^ + E - rA-v^mw dt 
•'^ i=i i=i 1=1 

<(« + l)e'(^lk(0)|Ui + ^|lE«;.(0)|U.^ <q|u(0)-i;(0)|U. (4.34) 
which together with the obvious inequality 

(n n 71 \ 

ii/^(E /i(E^^' - + E - /:(-f.)iiL^ 
i=l i=l j=l / 

finishes the proof of estimate (|4.28p . Passing now to the limit e ^ in that estimate (analo- 
gously to the proof of Theorem 14. ip . we sec that the limit semigroup S{t) generated by the ob- 
stacle problem p.ip satisfies indeed Assumption Jff with Bq := $. Theorem ()4.6p is proved. I 
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Corollary 4.7. Let the assumptions of Theorem \2.2\ hold and let the set Jlf he an n-dimensional 
simplex (|4.18p . Then, the semigroup S{t) associated with the obstacle problem (|l.ip possesses 
an exponential attractor M. in the sense of Definition \3.3[ Moreover, the global attractor A has 
finite fractal dimension. 

Indeed, this assertion is an immediate corollary of Theorems 13.41 and 14.61 
4.3 Some generalizations 

Wc now discuss the applications of our method to more general problems. We start with the 
obvious observation that all the above results remain valid if we replace the term Xu in the 
left-hand side of equation (jl.ip by any sufficiently regular interaction function (7(m, x). Namely, 
consider the problem 

dtu- Ar,u + diKiu) + g{x,u) 3 0, (4.35) 

where g G C(ri, C"'^(K", R")) is an arbitrary interaction function. Then, the following result 
holds. 

Theorem 4.8. Let Jlf he a convex bounded set o/ M" containing zero with a smooth boundary 
(or let K he an n-dimensional simplex (|4.18p ) and let g g C(0, C"'^(K", K")) he an arbitrary 
(not necessarily a gradient!) non-linear interaction function. Then, the solution semigroup S(t) 
associated with equation (|4.35p possesses an exponential attractor Ai in the sense of Definition 
\8.kA Moreover, the fractal dimension of the global attractor is finite. 

Indeed, since -J^ is bounded, the solution u is also automatically bounded in L°° (and 
the same will be true for the solutions of the approximate problems ()2.10p if e > is small 
enough no matter how the regular interaction function g looks like). So, the term g(u, x) can be 
treated as a perturbation and the proof of Theorem 14. 81 repeats word by word the given proof 
for the particular case g{x,u) := —Xu. 

Remark 4.9. The function g{x,u) may even depend explicitly on the gradient VxU, namely 
g = g{x,u,y xu)- However, in that case we already need to impose some restrictions on the 
growth of g with respect to (since the obstacle potential controls only the L°°-norm of a 
solution and the control of its W^'°°-norm should be then additionally obtained). In particular, 
if g does not grow with respect to V^jU, i.e., 

\gix,u,\7xu)\ < Qi\u\) 

for some monotone function Q independent of a; and VxU, the proof of Theorem 14. 81 still repeats 
word by word the case oi g{u) = —Xu. However, our conjecture here is that the result remains 
true under the standard sub-quadratic growth restriction 

\g{x,u,\/xu)\ < 0(|m|)(1 + |V,m|«) 

with q < 2. 

Remark 4.10. As we have already pointed out. our method of estimating the Lagrange multi- 
pliers is strongly based on the maximum principle for the leading linear part of equation (jl.ip . 
For this reason, we are unable in general to extend it to the case of non-scalar diffusion matrix. 
However, we point out that for some particular convexes Assumption ^ can be verified also 
for diffusion matrices (say diagonal). This is the case of the convex 

^:=[0,L]", L>Q. 
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In this case, Assumption can be easily verified (just multiplying the fcth equation of (|2.10|1 
by sgn(w^ — w^) even in the case of diagonal diffusion matrix 

dtu-al\^u + dlK{u)-\u3Q, (4.36) 

with a — diag(ai, • • • , afe), > for i = 1, • • • , 

One more generalization can be obtained replacing the Laplacian Aj;U in equations 
by the quasi-linear second order differential operator 

yl(u) := div(a(|up)V:cti), u = ("i, ■ ■ ■ , "rOi |up := H h 

with some natural assumptions on the scalar diffusion function a. Then, it is not difficult to 
verify that the proofs of Assumption given above remain true and the associated solution 
semigroup possesses an exponential attractor under the assumptions of Theorem 14.81 

5 Approximations of exponential attractors: the case of 
a simplex 

The aim of this section is to show that the exponential attractor M. of the singular problem 
(jl.ip can be approximated by the sequence of exponential attractors M.^ of the regular equa- 
tions (j2.10p . For simplicity, we consider below only the case in which X is an n-dimensional 
simplex (the case of an arbitrary bounded convex set with smooth boundary can be considered 
analogously). To be more precise, the main result of this section is the following Theorem. 

Theorem 5.1. Let the assumptions of Theorem \4.4\ hold. Then, for any sufficiently small £ > 0, 
there exists an exponential attractor for the semigroup Se(t) : (f>e — s- $j associated with the 
approximation problem ()2.10p (and the case e = corresponds to the semigroup S{t) ~ So{t) : 
$ — > $ associated with the limit singular problem (jl.ip / Moreover, the following conditions 
hold: 

• The exponential attractors are uniformly bounded in C'^~'^ for any i/ > 0. 

• Se(t)Mi; C Me o-nd the rate of convergence is uniform with respect to e 0: 

distL--iSeit)^e,Me) <Ce-"\ £>0, (5.1) 
where the positive constants C and a are independent of e. 

• The fractal dimension also remains bounded as £ 0; 

dim/(7We, i°°) < C < oo, (5.2) 

where C is independent of s. 

• The family is Holder continuous at e = Q: 

distfr(A<e,Xo) < C£^ (5.3) 

where the positive constants C and k are independent of e and dist*^™ stands for the 
symmetric Hausdorff distance between sets. 
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Proof. The proof of this Theorem is based on the foUowing abstract result taken from [14] 
(see also [TU]). 

Proposition 5.2. Let M' and ^\ he two Banach spaces such that J^i is compactly embedded 
in Jif and let, for any e > there exists a bounded closed set Bg C , a map : ^ B^ and 
two (nonlinear) projectors : B^ ^ Bq and : Bq ^ B^ such that the following properties 
hold: 

1) Uniform smoothing property: for every /ii,/i2 G B^, 

\\Sehi - < L\\hi - h2\\.^, (5.4) 

where L > is independent of e > and h\,h2 G Bg. 

2) The maps and §o o,'"'^ close in the following sense: 

||§,oQ,;j-§o/i||,_^ < Ce, V/ieBo, 
\\%oTl,h-^,h\\,yf <Ce, V/ieB^, 

where the constant C is independent of e and h. 

Then, the discrete semigroups Se(n) generated by the maps possess a uniform family 
of exponential attractors C B^ such that the following properties hold: 

1) Uniform rate of attraction: 

dist^(§e(n)Be,7We) < Ce""" (5.6) 

for some positive C and a independent of e. 

2) Uniform bounds for the fractal dimension: 

dimf{M„Jf) <C, (5.7) 

where C is independent of e. 

3) Holder continuity at e = 0: 

distJ"(7W„7Wo) < Ce" (5.8) 
for some positive C and k independent of s. 

We can now prove Theorem 15.11 by combining (as we did in the proof of Theorem 13. 4p 
the Proposition above with the £ trajectory approach. Namely, we define the spaces and 
Jffi by (|3.14p and introduce the trajectory phase spaces 

B. := {ueJ^i, u{t) = Se{t)uQ, t e [0, 1], uo e $£} (5.9) 

and, using the lifting solution operator 

Te:$e^B„ {T.uom -.^ S,{t)uo, 

we lift the solution semigroup S'e(t) to the trajectory phase space B^: 

§e(t) := T, o S,{t) o T-\ S,it) : Be ^ Be 

and set, finally, := Se(l). 

Let us verify the assumptions of the abstract proposition for the maps thus defined. 
To this end, we note that, thanks to Theorem 14.61 we have the uniform estimate (|4.28|1 for 
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the approximations to the Lagrange multiphers — /i„2 associated with two trajectories 
ul Se{t)uQ, i ~ 1,2. Thereofere, arguing exactly as in the proof of Lemma 13.71 we conclude 
that 

- "?llL2(i,2;ffi) + \\dtul - 9tUe||Li(l,2;ff-») < - U^||l2(o,1;L2) (5.10) 

for some positive L which is independent of s and ul, i ~ 1,2. Thus, assumption (|5.4p is 
verified. 

In order to verify (jS.Sp . we first fix an arbitrary interior point wq G of Kq :~ K 
and introduce the linear contraction map Eg, : K (where is defined by (|4.20|) ) by the 

following expression: 

1 1 re 

Ee(.u)-~— {u-Wq) + Wq^— U+— Wo, (5.11) 

1 + re 1 + re 1 + re 

where r is a sufficiently large (but independent of e) positive constant. Indeed, it is not difficult 
to show using the explicit expression for the set Kg, and for the function F^, that E^{K^) C A'o 
if r > is large enough. Finally, we define the map : ^ ^ as follows 

Ee{u){x) -.^ (pe{x)Es{u{x)), 

where the cut-off function (pe(,x) € C^(R^) is such that < (p{x) < 1, (pe{x) = if x G dil and 
ipe{x) = if dist(.T,9ri) > e. Then, obviously, 

\\il-E,)v\\L2<Ce\\v\\L2, yve^e (5.12) 

and for some positive C independent of e and v. 
We now need the following crucial lemma. 

Lemma 5.3. Let the assumptions of Theorem \5.1\ and let E^ he defiend by (|5.1ip . Then, for 
every u^ G and every uq G Kq, the following estimate holds: 

\\S,{t)u^, - Soit)uo\\L^ < C(£i/2 + _ u^||L2)e^*, (5.13) 

where the positive constants C and K are independent of e, G and uq £ $. 

Proof. Let u{t) ■= So{t)uo, Ue{t) Ss{t)u° and Ve{t) := Ee{ue{t)). Then, since v^{t,x) £ 
K for alH > 0, it is an admissible test function for the variational inequality (|2.2p . Therefore, 

{dtu{t), u{t) - v,{t)) + (V^S(t), V:,ii(t) - V:,ue(t)) < X{u{t),u{t) - v,{t)) (5.14) 

for almost all t>Q. Introducing the function We{t) ■— u^{t) — v^{t) = (1 — E^)ug{t) and ()5.12p . 
we have 

{dtu{t), u{t) - u,{t)) + (V^?2(t), V:rw(t) - V:,ue(i)) - A(?2(t), u{t) - u,{t)) = 

= {dtu{t), u{t) ~ v,{t)) + (V,il(t), V,u{t) - ^^xVeit)) - X{u{t),u{t) - v,{t))+ 
+ {dtu{t), u,{t) - v,{t)) - {/\,u{t),u,{t) - v,{t)) - \(u(t), u,(t) ~ v,(t)) < 

<-{hu„{t),We{t))<Ce, (5.15) 

where we have implicitly used that the Lagrange multiplier h.u{t) associated with the solution 
u{t) of the singular problem (|2.20p is uniformly bounded (see Proposition 12. 5p . Moreover, 
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multiplying equation (j2.10l) by u^{t) — u{t), integrating over x € fl and using the monotonicity 
of and the fact that F^{u{t)) = 0, we arrive at 

{dtu,{t),ue{t) - u{t)) + {V,^u,{t),V^Ueit) - V^u{t)) - X{u,{t),u,{t) - u{t)) < (5.16) 

which holds for almost alH > 0. Taking a sum of (jS.lSp and (|5.16p . we finally have 

~\K{t) - uimh + \\^.u,{t) - \/Mmh - MKit) - u{t)\\l, < Ce (5.17) 

which, together with the Gronwall inequality, gives ()5.13p and finishes the proof of the Lemma. 

I 

It is now not difficult to finish the proof of the Theorem. To this end, we set Ilg : ^ Bo and 
: Bo ^ Be by 

m := To o o T7\ :=T,oT(7i, 

where Tq is the lifting operator related to SQ{t). Then, it is not difficult to see using Lemma 
15.31 and estimate (|5.12p that the projectors 11^ and Qe satisfy estimates (|5.5p . Thus, all of the 
assumptions of the abstract Proposition 15.21 are verified and, consequently, the discrete semi- 
groups §e(n) acting on the trajectory phase spaces possess the uniform family of exponential 
attractors which satisfies conditions (|5.6p , ()5.7p and ()5.8p of Proposition 15.21 

The rest of the proof can be made exactly as in Theorem 13.41 Indeed, arguing as in 
Lemma ITU we see that 

\\ul{l) - ul{l)\\l. < (2A + 1) /' \\ul{t) - ^.2(^)111^ (5.18) 

Jo 

for any two solutions u\ and of the approximation problems (|2.10p . In addition, multiplying 
(|5.17p by t and integrating over t <S [0, 1], we see that 

||u,(l)-i2(l)||2, < (2A + 1) / \\u,{t)-u{t)\\l,dt. (5.19) 

Thus, projecting the constructed exponential attractors back to the physical phase spaces 
$e by 

we obtain the uniform family of exponential attractors for the discrete semigroups (n) : — > 
$e which satisfy (|5.ip . (|5.2p and (|5.3p in a weaker space instead of In addition, we have 
also that 

and, therefore, due to estimate (|4.23p and to Corollary 12.71 these attractors are uniformly 
bounded in C^~'^{n) for all u > and the maps (i,Mo) — > Se(t)uo are uniformly Lipschitz 
continuous on [0, 1] x M^: (with the Lipschitz constant independent of e as well). Thus, the 
standard formula 

Me -.^ Ute[OA]Seit)Mf 

gives the desired uniform family of exponential attractors for the semigroups S^it) : $£ ^ $£ 
with continuous time. This family of exponential attractors satisfies ()5.ip . (|5.2p and (j5.3p in a 
weaker topology of L^(0) instead of L°°(0) but, using the fact that Ai^ are uniformly boimded 
in C^~'^(r2) together with the interpolation inequahty ()3.19p . we see that ()5.ip . (|5.2p and ()5.3p 
hold in the initial topology of L°° as well. Theorem 15. II is thus completely proved. | 
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Remark 5.4. It is not difficult to show that, under the assumptions of Theorem 15.11 the 
solution semigroups Se{t) are actually defined not only on the space $£, but in much lager 
phase spaces, namely, on L^. Moreover, arguing analogously to the proof of Theorem 14.41 one 
can see that the set $e is an absorbing set for the semigroup S'e(t). This means that, for every 
bounded B C L'^{fl) there exists T = T{\\B\\) (independent of e) such that 

^,(t)Bc$„ t>T. 

Thus, the constructed exponential attractors attract (uniformly with respect to e) the 
bounded sets of L^{U) as well. 

Finally, one may even extend the solution semigroup S{t) associated with the obstacle 
equation (jl.ip on the whole space as well. To this end, we just need to pass to the limit 
e ^ in the solutions of the approximate problems (|2.10p . It is not difficult to show that the 
associated solution u{t) := S{t)uQ will belong to $ for every t > and, in the case uq ^ $, it 
will have a jump at t = and the i^-limit 

Uq := lim u{t) e <I> 

will exist. Thus, we factually have 

u{t) = S{t)uo - 5(t)n$(uo), 

where 11$ : L'^{fl) ^ $ is a non-linear "projector" to the set (f>. 

However, in contrast to the semigroup S(t) : $ <f> which is uniquely defined by 
the singular equation p.ip . the above non-linear "projector" 11$ essentially depends on the 
particular choice of the approximation of the obstacle potential and, therefore, is not canonically 
defined by the problem (jl.ip itself. For this reason, we have preferred not to use such projectors 
in our paper and to define the solution u{t) for the initial data u(0) € only. 
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